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ABSTRACT 

We use numerical simulations to investigate the orientation of the angular momentum axis of disk 
galaxies relative to their surrounding large scale structure. We find that this is closely related to the 
spatial configuration at turnaround of the material destined to form the galaxy, which is often part of 
a coherent two-dimensional slab criss-crossed by filaments. The rotation axis is found to align very well 
with the intermediate principal axis of the inertia momentum tensor at this time. This orientation is 
approximately preserved during the ensuing collapse, so that the rotation axis of the resulting disk ends 
up lying on the plane traced by the protogalactic material at turnaround. This suggests a tendency for 
disks to align themselves so that their rotation axis is perpendicular to the minor axis of the structure 
defined by surrounding matter. One example of this trend is provided by our own Galaxy, where the 
Galactic plane is almost at right angles with the supergalactic plane (SGP) drawn by nearby galaxies; 
indeed, the SGP latitude of the North Galactic Pole is just 6 degrees. We have searched for a similar 
signature in catalogs of nearby disk galaxies, and find a significant excess of edge-on spirals (for which the 
orientation of the disk rotation axis may be determined unambiguously) highly inclined relative to the 
SGP. This result supports the view that disk galaxies acquire their angular momentum as a consequence 
of early tidal torques acting during the expansion phase of the protogalactic material. 

Subject headings: Galaxy: disk, structure, formation 



at early times (since f2 sa 1 then) and effectively ends at 
turnaround, when I is maximal. In general, then, the di- 
rection of the angular momentum will be determined by 
the shape of the protogalactic material at turnaround, and 
is expected to align with the intermediate axis of inertia if 
T and I are uncorrelated. 

Many of these expectations have been confirmed by a 
number of studies of the angular momentum properties 
of dark matter halos formed in cosmological simulations 
(see, e.g., Barnes & Efstathiou 1987, as well as the more 
recent work of Porciani, Dekel & Hoffman 2002a, b, and 
references therein). The latter authors point out that T 
and I are, in fact, highly correlated, and conclude that, 
relative to the principal axes of T, the most significant 
trend is for the spin to be perpendicular to the direction 
of maximum compression (i.e., that of the minimum eigen- 
vector of T). Because T and I are correlated, the direction 
of maximum compresion roughly coincides with the direc- 
tion of the minor axis of inertia at the time of turnaround. 
A solid prediction of TTT is, therefore, that galaxy spins 
should be nearly perpendicular to the minor axis of the 
collapsing material at late times. 

The main shortcoming of this work is that it applies 
mainly to dark halos, rather than to the luminous (ob- 
servable) component of galaxies. Simulations that include 
the presence of a dissipative gaseous component have high- 
lighted the possibility that large losses of angular momen- 
tum may accompany the collapse of the baryonic com- 
ponent (Navarro & Benz 1991, Navarro & White 1994, 
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1. INTRODUCTION 

In hierarchical models of galaxy formation, the origin 
of galactic angular momentum is ascribed to tidal torques 
operating early on the material destined to form a galaxy. 
Over the years, starting from the ideas of Stromberg (1934) 
and Hoyle (1949), and quantified by the work of Peebles 
(1969), Doroshkevich (1970), and White (1984), a number 
of important properties inherent to the acquisition of an- 
gular momentum through tidal torques (hereafter referred 
to, for short, as "tidal torque theory", or TTT) have been 
identified. 

The efficiency of the torquing is low, typically endow- 
ing galactic material with a rather small amount of net 
rotation. This implies that large collapse factors within 
massive dark matter halos are needed to explain the 
centrifugally-supported nature of galactic disks (Fall & Ef- 
stathiou 1980). To leading order, galaxy spins result from 
the misalignment between the principal axes of the iner- 
tia momentum tensor (Iy) of the material being torqued 
and of the "shear" tensor (Ty = — d 2 4>/ dxidxj) gener- 
ated by external material. In the Cartesian principal axis 
frame of the protogalaxy, the leading term of the torque is 
given by t% — dLi/dt ~ Tjk{Ijj — Ikk) (here i, j, and k are 
cyclic permutations of 1 to 3, and Li are the cartesian com- 
ponents of the angular momentum). The inertia term is 
maximal along the direction that maximizes the difference 
between Ijj and i-e., the intermediate axis of inertia. 
Angular momentum growth is typically linear with time 
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All baryons (r<2.5 comoving Mpc) at z ~0 



Galaxy baryons aL turnaround (z ~3) 
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Fig. 1. — (a-left) Spatial distribution of all baryons in the resimulated high-resolution region at late times. Three orthogonal projections 
are shown, aligned with the principal axis of inertia of the system at that time. Units are in physical kpc. Particles in yellow are stars; others 
are gas. Labels in each panel indicate the aspect ratio of the system. Arrows indicate the relative importance of the cartesian components of 
the angular momentum. Bottom-right panel is a zoom of the Y-Z projection that shows the gaseous disk almost face-on, implying that its 
rotation axis lies on the plane of the surrounding structure, (b-right) Same as (a) but for the baryons that will collapse to form the galaxy 
at 2 = 0, shown at turnaround in its own principal axis frame. The bottom-right panel shows a zoomed-in (X',Z') projection of the same 
baryons at 2 ~ 0, which have collapsed to form a disk whose rotation axis remains closely aligned to the intermediate (V) axis at turnaround. 



Navarro & Steinmetz 1997). The rather indirect mapping 
between the angular momentum of dark halos and that of 
their baryonic components makes it difficult to assess the 
success of TTT in accounting for the spin of spiral galax- 
ies. The discussion of the previous paragraphs, however, 
suggests that TTT may be tested by looking for resid- 
ual correlations originating from the coupling between the 
shear and inertia momentum tensors that dominates the 
angular momentum growth in protogalaxies. One example 
of this is the recent work of Lee & Pen (2000, 2001, 2002, 
see also Lee, Pen & Seljak 2000), who attempt to recon- 
struct the shear field from the observed galaxy spin fields 
and to uncover correlations between the spin field and the 
large-scale distribution of matter. 

We investigate here a related signature expected from 
TTT; i.e., the correlation between the "shape" of proto- 
galactic matter distribution at turnaround and the direc- 
tion of the resulting angular momentum. As explained 
above, TTT suggests a trend for spin axes to be perpen- 
dicular to the minor axis of the sheet, so that disk galaxies 
should be highly inclined relative to the plane defined by 
their surrounding structure. 

In this Letter we explore whether disk galaxies assem- 
bled hierarchically in ACDM numerical simulations indeed 
follow these expectations. In addition, we search for a 
statistical signature of the predicted trend in catalogs of 
nearby galaxies. 

2. NUMERICAL EXPERIMENTS AND RESULTS 

We have analyzed four N-body/gasdynamical simula- 
tions of the formation of disk galaxies in the "concor- 
dance" ACDM cosmogony (fi = 0.3, n A = 0.7, h = 0.65, 



Q b = 0.019 h- 2 , a s = 0.9). These simulations are simi- 
lar to those described in detail by Abadi et al (2003a, b) 
and Meza et al (2003), where we refer the reader for full 
details. The simulations follow self-consistently the evolu- 
tion of a small region surrounding a target galaxy, excised 
from a large periodic box and resimulated at higher reso- 
lution preserving the tidal fields from the whole box. For 
the purposes of the discussion here all four simulations 
give consistent results. We have chosen to illustrate these 
using the simulation presented by Abadi et al (2003a, b), 
but will use the others to verify the general applicability 
of our results. 

Figure la shows, at late times (z ~ 0), the baryons 
within ~ 2.5 Mpc of the target galaxy. The three dif- 
ferent projections in the figure have been chosen so as to 
coincide with the principal axes of the inertia momentum 
tensor. This figure shows the high degree of anisotropy 
that characterizes the non-linear evolution of structures in 
a ACDM universe. What started off as a roughly spher- 
ical region ~ 5 comoving Mpc across develops into a co- 
herent triaxial structure that surrounds the target galaxy 
(located at the center of the panels in Figure la). The mi- 
nor axis of the structure, in particular, is very stable across 
the structure, which may thus be described as a "sheet" 
criss-crossed by filaments. The ratios of the inertia tensor 
eigenvectors are given in Figure la, and are not dissimilar 
to the mean of all four systems: (I22/I11) — 0.419 (with 
dispersion a = 0.360), (hs/hi) = 0.174 (a = 0.157) and 
(I33/I22) = 0.408 (a = 0.059). 

The arrows in the panels of Figure la indicate the rel- 
ative magnitude of the cartesian components of the angu- 
lar momentum of the sheet, and illustrate the TTT pre- 
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Fig. 2. — (a- left) Histogram of supergalactic position angles of edge-on (b/a < 0.175) spirals in the Principal Galaxy Catalog (PGC, Paturel 
et al 1997) selected within various recession velocity limits, as labelled. A position angle of 0° means that the galaxy's plane is perpendicular 
to the SGP; 90° means that it is parallel to the SGP. There is a well-defined excess of spirals perpendicular to the SGP in the vicinity of the 
Milky Way. (b-right) Aitoff equatorial projection of all spirals within 1200 km s — 1 in the PGC. The U-shaped thick curve is the Galactic 
plane (GP); the S-shaped curve is the SGP. The projected major axis of edge-on spirals is shown for galaxies with position angles smaller 
than 35° (filled circles) and greater than 55° (open circles). 



dieted trend (see §1) for spins to be approximately per- 
pendicular to the direction of maximum compression or, 
roughly, to that of the minor axis (Z) during the non- 
linear stages of the evolution. Indeed, the whole sheet of 
material is slowly spinning around the intermediate prin- 
cipal axis of inertia (i.e., \Ly\ 3> \Lx\, \Lz\)- This is a 
generic result of our four simulations; quantitatively, we 
find that the angular momentum of the high-resolution re- 
gion surrounding the target galaxies to be oriented such 
that (\L X \/\L\) = 0.255 {a = 0.104), (\L Y \/\L\) = 0.860 
(cr = 0.103), and (\L Z \/\L\) = 0.351 (a = 0.275). 

This two-dimensional structure is present since early 
times and, at z ~ 3, its minor axis is approximately paral- 
lel to that of the material that will collapse to form the tar- 
get galaxy at z = (i.e., the "protogalaxy" , see Figure lb, 
as well as Sugerman, Summers & Kamionkowski 2000 for 
a similar analysis). The shape of the protogalactic mate- 
rial at turnaround is indicated in the the panels of Figure 
lb; for the four simulations we find (I22/I11) = 0.652 (with 
dispersion a = 0.0.247), (hs/hi) = 0.265 (cr = 0.095) and 
(I33/I22) = 0.429 (cr = 0.132). The angular momentum of 
the protogalaxy at z ~ 3 shows similar alignment proper- 
ties to that of its surrounding material at z ~ 0; in par- 
ticular, the close alignment with the intermediate inertia 
axis at turnaround is well defined in all four simulations: 
(\L X ,\/\L\) = 0.055 (a = 0.031), (\L Y ,\/\L\) = 0.832 
(cr = 0.260), and (\L Z ,\/\L\) = 0.421 (cr = 0.319). 

The angular momentum direction of the baryons at 
turnaround is approximately preserved by the disk com- 
ponent at z = and, as a result, the plane of the disk 
is highly inclined relative to the plane defined by its sur- 

5 We define as turnaround the time when the moment of inertia of 



rounding large-scale structure. This is illustrated in the 
bottom right panels of Figures la and lb, where the disk 
at z — is seen approximately face-on in projections that 
contain cither the minor axis of the surrounding struc- 
ture at z ~ or that of the protogalaxy at turnaround. 
This indicates that the shape of the protogalactic mate- 
rial during the expansion phase effectively determines the 
orientation of the spin axis of the resulting galaxy. Since 
these anisotropies correlate across a wide range of scales 
in ACDM, a signature of this process might be preserved 
as a residual alignment between the disk and its surround- 
ings; the rotation axis of disk galaxies may be expected to 
lie on the plane defined by its surrounding structure. We 
explore next whether such alignments are indeed present 
in samples of disk galaxies in the nearby Universe. 

3. ORIENTATION OF NEARBY DISK GALAXIES 

As recognized by de Vaucouleurs (1953) in the Shapley- 
Ames catalog, galaxies in the vicinity of the Milky Way 
are arranged in a two-dimensional slab usually referred 
to as the "super galactic plane" (SGP). This structure 
stands out clearly in whole-sky catalogs of nearby galax- 
ies; extends out to several tens of Mpc from the Galaxy; 
and includes many of the major features of the local large 
scale distribution of matter, such as the Virgo cluster, the 
Perseus-Pisces supercluster, and the Great Attractor (La- 
hav et al 2000). 

The plane of the Galaxy is highly inclined relative to the 
SGP; indeed, the two planes are approximately perpendic- 
ular to each other, as indicated by the low supergalactic 
latitude of the North Galactic Pole (~ 6°). This situation 
is similar to that of the simulated disk galaxy in relation 

protogalaxy attains its maximum. 
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to its surrounding structure, as discussed in §2. It is there- 
fore tempting to regard the rather peculiar orientation of 
the Galactic plane relative to the SGP as a result of early 
torques acting during the protogalactic stage. 

If this interpretation is correct, we would expect an ex- 
cess of nearby galaxies whose rotation axes are approxi- 
mately perpendicular to the normal to the SGP. The ex- 
istence of such alignment has proven controversial; prior 
work has led to claims of statistical evidence for such "anti- 
alignment" of galactic planes relative to the SGP (see, e.g., 
Flin & Godlowski 1989 and references therein), but also 
to suggestions that there is no convincing evidence for de- 
viations from a random distribution of orientations (see, 
e.g., Dekel 1985) 

A number of factors complicate this analysis and may 
explain this disagreement. For example, the statistics of 
the position angles of galaxy samples that do not discrim- 
inate properly between elliptical and disk galaxies may 
be inconclusive, as the projected major axis of early type 
galaxies — in contrast to those of disks — may bear little 
relation to the angular momentum of the galaxy. In ad- 
dition, determining the direction of the rotation axis of a 
spiral galaxy requires knowledge not only of the position 
angle and inclination (readily available in most catalogs), 
but also of how the galaxy is inclined on the sky, i.e., which 
side of the galaxy is closer to the observer. This is known 
for a few well-studied spirals, such as M31 6 , but it is 
generally unavailable for most galaxies. The uncertainty 
introduced by this ambiguity depends on inclination: it is 
maximal for galaxies inclined by 45 degrees, but is neg- 
ligible for galaxies seen either face-on or edge-on. In this 
Letter we restrict our analysis to edge-on spirals and search 
for a possible correlation between the orientation of galaxy 
disks and the SGP. We explore ways of circumventing the 
aforementioned ambiguity and extend this analysis to all 
nearby spirals in a forthcoming paper. 

Figure 2a shows the distribution of position angles (in 
the supergalactic reference frame) of all edge-on nearby 
spiral galaxies in the Principal Galaxy Catalog (PGC) 
with recession velocities within various limits. Since we 
are only interested in the direction of the rotation axis rel- 
ative to the SGP, we have reduced all position angles to 
between 0° and 90°. As shown in Figure 2a, there is a 
clear excess of nearby edge-on galaxies highly inclined rel- 
ative to the SGP. The significance of the excess decreases 
the larger the volume considered around the Milky Way. 
Within 1200 km s _1 there are 196 spirals in the PGC with 
measured inclinations and position angles. Of these ~ 30 
are edge-on. Figure 2a shows that roughly three times as 
many edge-on galaxies with position angles between 0° and 
20° as in the range (70°, 90°); a simple KS test shows that 
this excess is significant to the 92% level. The magnitude 



of the excess (and its significance) decreases to ~ 40% for 
spirals with recession velocities under 2400 km s _1 , and 
is essentially negligible when galaxies within 5000 km s _1 
(or greater) are considered. 

This result supports the tendency of spirals to be highly 
inclined relative to the SGP claimed by Flin & Godlowski 
(1989). Our study provides a compelling physical inter- 
pretation for the origin of such alignment and supports 
the view that the angular momentum of spiral galaxies 
originates in tidal torques tightly coupled to anisotropies 
during the expansion, turnaround, and collapse of proto- 
galactic material. 

4. SUMMARY AND DISCUSSION 

We have used cosmological N-body/gasdynamical sim- 
ulations to study the orientation of disk galaxies relative 
to their surrounding structure. We find that the angu- 
lar momentum axis of simulated galaxies is perpendic- 
ular to the minor axis (and very well aligned with the 
intermediate axis) of the material destined to form the 
galaxy at turnaround, which is usually arranged on a two- 
dimensional slab criss-crossed by filaments. These sheet- 
like structures may extend out to several Mpc, in a way 
reminiscent of the "supergalactic plane" that surrounds 
the Milky Way. This coherence over a wide range of 
scales in the anisotropies present during the expansion 
and turnaround of protogalactic material generally results 
in galaxy disks highly inclined relative to the large-scale 
two-dimensional structure where they are embedded. 

This provides a natural explanation for the high inclina- 
tion of the Milky Way relative to the supergalactic plane, 
as well as for the excess of nearby edge-on spirals whose 
rotation axes lie approximately on the supergalactic plane. 
It may also offer a physical interpretation for the tendency 
of satellite galaxies to align along the minor axis of bright 
spirals (the "Holmberg effect"). Our interpretation offers 
as well a number of "natural" predictions that may be 
used to falsify it. In particular, we expect an excess of 
highly inclined spirals in any two-dimensional large-scale 
distribution of galaxies, a result that may have already 
been detected in the Perseus supercluster (Flin 1988), and 
that could be verified in large galaxy surveys such as the 
2dfGRS and the SDSS. The detection of these and other 
non-trivial correlations between the spin and matter fields 
would serve to establish beyond doubt the validity of tidal 
torque theory as the origin of the angular momentum of 
spiral galaxies. 

We thank D.Garcia Lambas and L. Sales for useful dis- 
cussions, and the referee, Avishai Dekel, for a constructive 
report. 
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